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^^ ' We report the first observation of exclusive decays of the type B — > D*NNX, 

•^ ■ where A^ is a nucleon. Using a sample of 9.7 x 10® BB pairs collected with the 

CLEO detector operating at the Cornell Electron Storage Ring, we measure 
the branching fractions B (S°^ D*" p p 7r+) = {d.BtH ± 1.0) x 10""^ and 
B{ B^^ D*- pn) = (14.5t|o ± 2.7) x 10^"^. Antineutrons are identified by 
^^ ' their annihilation in the Csl electromagnetic calorimeter. 
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A unique feature of the B meson system is that the large mass of the b-quark allows for 
many of the weak decays of the B meson to include the creation of a baryon-antibaryon pair. 
In the simplest picture, baryons are expected to be produced in decays of the type B -^ 
AcpX, and it has been only decays of this type which have been exclusively reconstructed 
to date [0. However, one can combine the recently measured value 13{A^ — > pK^n^) = 
(5.0±0.5±1.2)% with estimates of the product branching fraction B{B -^ A^X) xB{Ac -^ 
pK'n+) of (1.81 ± 0.22 ± 0.24) x lO^^ [|] to determine that B -^ A^NX modes, where A^ 
is a proton or a neutron, account for only about half of the total B -^ Baryons rate. 
Dunietz 0] has suggested that modes of the type B —>■ DNNX, in which D represents any 
charmed meson, are likely to be sizeable. B -^ DNNX final states can arise from either 
the hadronization of the W boson into a baryon-antibaryon pair, or the production of a 
highly excited charmed baryon which decays strongly into a baryon plus a charmed meson. 
CLEO has previously reported an inclusive upper limit for B{B -^ DNNX) of < 4.8 % at 
90% confidence level [^. We report the first observation of decays of this type, and present 
measurements of the branching fractions B{B'^^ D*^ p p n^) and B{B'^^ D*^ p n). The 
charge conjugate process is implied in the reconstruction of B^—>- D*^ p p n^. However, 
in the reconstruction of B^^ D*^ p n only the mode with the antineutron is used in our 
measurement because neutrons do not have a distinctive annihilation signature in the CLEO 
detector. 

The data were taken with the CLEO II detector at the Cornell Electron Storage Ring 
(CESR). The sample we use corresponds to an integrated luminosity of 9.1 fb~^ from data 
taken on the T(4S') resonance, corresponding to 9.7 x lO^BB pairs, and 4.5 fb~^ in the 
continuum at energies just below the T(4S'). We assume that 50% of the BB pairs consist of 
B^B^, and that there are equal numbers of i?° and B^ mesons. Charged particle trajectories 
are measured in a cylindrical drift chamber operating in a 1.5 T magnetic field. Photons and 
antineutrons are detected using an calorimeter consisting of 7800 Csl crystals with excellent 
resolution in position and electromagnetic shower energy. Simulated events were generated 
with a GEANT-based Monte Carlo program |^. Sixty percent of the data were taken in the 
CLEO II. V configuration [|. 

Charged particle identification is accomplished by combining the specific ionization 
(dE/dx) measurements from the drift chamber with time-of-flight (TOP) scintillation 
counter measurements. We reconstruct the decay mode D*^ -^ D^tt^, with D^—*K~^7r~ , 
D'^ ^ K~^ n~ TT^ , and D^—>-K~^n~'K~^'K~. Pairs of calorimeter showers with photon-like lateral 
shower shapes and invariant mass within 2.5 standard deviations of M(7r°) are considered as 
7r° candidates. We select D*~ candidates using 95% efficient cuts around the central values 
of the Mdo and (Md* —Mdo) distributions. In the few cases where there is more than one D* 
candidate in an event, the candidate which is closest to M^io = 1.8646 GeV and {Md* —Mdo) 
= 0.1454 GeV g is chosen. 

To reconstruct B^ candidates for B^ -^ D*^ p p n^ , we calculate the beam constrained 



mass, M(B) = y -Ebeam^ — p(-B)^, where -Ebeam is the beam energy, and p{B) is the B^ 
candidate three-momentum magnitude. We also use the energy difference between the beam 
energy and the energy of the reconstructed B candidate: AE = £'beam — E{B). Using a 
Monte Carlo simulation program, we find the detector resolution for the AE distribution for 
each D^ mode, and require that the measured AE is within 3 standard deviations of zero. 
The M{B) resolution is dominated by the beam energy spread and is consistent with the 



Monte Carlo simulation for B'^ -^ D*^ p p tt^, which predicts it to be Gaussian with a width 
of cr = 2.7 MeV. The M{B) distribution is fitted to a Gaussian function with the predicted 
width, and a polynomial background function with suppression at the -Ebeam threshold. The 
fitted signal yield is 32.3tg;Q events, where the errors are statistical only. If a is allowed to 
float, a value of a = 2.1 ± 0.4 MeV is obtained, consistent with the Monte Carlo prediction. 
Figure |l] shows AE vs M{B) for B^ -^ D*^ p p tt^, and Figure |](a) shows M{B) for B^ -^ 
D*~ p p 7r+. 
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FIG. 1. AE vs M{B) distribution for S°- 
the signal region. 
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FIG. 2. M{B) distribution for (a) B°^ D*' p p Tr+ and (b) B^^ D*' p n. Each plot is fitted 
to the sum of a background function and a fixed-width Gaussian signal shape. 

To find the decay B^—^ D*~ p n requires the detection of an antineutron. This is accom- 
plished by identifying annihilation showers in the Csl calorimeter. Antiprotons and antineu- 
trons annihilate with nucleons in the calorimeter. These annihilations result in showers with 



different characteristics than those from photons or charged particles. We use antiproton 
annihilation showers to define the antineutron selection criteria since we are unable to iso- 
late a sample of antineutrons in data. Our antiproton sample consists of 1.6 x 10^ p's from 
reconstructed A -^ p7r~^ decays, in which the daughter antiprotons are selected by dE/dx 
and TOF response. The isolation of this sample is independent of calorimeter response and 
therefore allows us to evaluate our shower-based selection criteria. We can expect these same 
criteria to be similarly effective for antineutrons. 

We find that antinucleons typically deposit a substantial amount of energy in a laterally 
broad shower, which has energy i^main, and a lesser amount of energy in adjacent satellite 
showers, whose energy is added to -Emain to define -Egroup- The selection of showers with -Emain 
> 500 MeV and -Egroup > 800 MeV is useful in suppressing backgrounds, while retaining many 
antinucleon annihilation showers. The main shower must have polar angle 6 with respect 
to the incoming positron direction of 45° < 9 < 135° to ensure that it is located in that 
part of the calorimeter that has the best resolution, must have a lateral shape broader 
than typical photons, and must not match the projection of any charged track trajectory. 
Baryon number conservation serves as an added background suppressant: once we require 
that there be a proton, the likelihood of an annihilation-like shower in the event to be an 
antineutron increases significantly. The reconstruction efficiency as a function of momentum 
for antineutrons in a Monte Carlo sample of B^^ D*~ p n events in which the B^ selection 
criteria have been applied is shown in Figure ^ 
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FIG. 3. Antineutron reconstruction efficiency derived from simulated B^ — > D* pn events. 



The measured shower energy for an antineutron, E'group, does not give an accurate mea- 
surement of the total energy of the antineutron. We therefore assign the energy of the 
antineutron candidate using E{n) = .Ebeam — E{D*^) — E{p). We then use this energy, 
together with the position of the antineutron shower, to calculate the momentum of the an- 
tineutron candidate. This momentum is added to the momenta of the D*~ and p candidates 
to give the momentum of the B^ meson, p{B). The B^ candidate mass is then calculated 
with M{B) = ^El^^-p{Bf. 

The mass resolution of the reconstructed B^ meson is estimated from a Monte Carlo 
simulation to be Gaussian with a width a = 3.1 MeV, demonstrating that the lack of a 
direct measurement of the antineutron energy does not seriously degrade the M{B) resolution 
relative to B^ — > D*^ p p n^. However, for this mode we cannot use the requirement that 



AE is consistent with zero, as we do not have two independent measures of the energy of 
the 5°. 

In reconstructing 5*^— > D*^ p n we could also be reconstructing i?° — > D*^ D^ with Df 
-^ p n. This latter decay has not been observed but could occur via cs annihilation. We 
might also be including events that are S° -^ D*^ D*'^ with £>*+ -^ Dfj or £>*+ -^ D^t^^ 
and D"^ -^ p n. These events will populate the M{B) signal region but with a broader 
signal peak due to the missing soft photon or 7r°. We eliminate these two types of decay by 
rejecting events with 1.91 < M{pn) (GeV) < 2.04 GeV, for a loss of only 9% in the relative 
reconstruction efficiency. 

The final data M{B) distribution, shown in Fig. 0(b), is fit with a Monte Carlo predicted 
Gaussian width of cr = 3.1 MeV and gives a signal yield of 24.0l|;o events. If a is allowed to 
fioat, a value of 2.6 ± 0.4 MeV is found, consistent with the Monte Carlo expectation. 

We use a Monte Carlo simulation to calculate detection efficiencies to be 0.52% for i?°— >• 
D*^ p p TT"*", and 0.34% for B^^ D*~ p n, where these numbers include all the relevant 
branching fractions of the daughters. The Monte Carlo generation of the decays assume no 
resonant sub-structure. No evidence of sub-structure has been found in the signal candidates. 
However, we will allow for possibility of sub-structure, which would alter the efficiency, in 
our estimation of the systematic uncertainties. 

Whereas the Monte Carlo simulation has been well tested for charged particles and pho- 
tons, this is the first measurement that has explicitly needed the efficiency for antinucleon 
annihilations in the calorimeter. We find a discrepancy for antiproton annihilation showers 
between the Monte Carlo and data. The reconstruction efficiency for Monte Carlo and data 
for antiprotons, as determined from the aforementioned A sample, is shown in Fig. ^a). We 
find that the antiproton reconstruction efficiency is overestimated by the GEANT simulation 
prediction. We assume that the Monte Carlo simulation for the antineutron efficiency is sim- 
ilarly overestimated and therefore must be corrected. Weighting the efficiency correction for 
antiprotons by the expected antineutron momentum spectrum results in an antineutron se- 
lection efficiency which is 21% lower than the GEANT Monte Carlo simulation. In Fig. ^(b) 
we show the antineutron momentum spectrum generated in Monte Carlo for B —^ D*~pn 
decays. 

We search for, but do not find, resonant substructure contributions to B^ -^ D*^pp'K^ 
and B^ -^ D*^pn decays. Examples of possible substructure arise from a heavy charmed 
baryon decaying strongly to {p D*~) for B^^ D*~ p p -k^ and (n D*~) for B^^f D*~ p n, 
and a resonance of the virtual W decaying to ppn^ . We also study the effect on the Monte 
Carlo reconstruction efficiency of a two-body decay into D*~X, and find no evidence from 
the kinematic distribution of the daughter particles for decays of this type. We also find 
insignificant background from B decays with a charmed baryon final state or from other 
B -^ D*-X modes. 

We have considered many sources of systematic uncertainty of these branching fraction 
measurements. Systematic uncertainties shared by both modes are statistical uncertainty of 
D^ branching fractions (0.6%) and D*^ branching fraction (1.4%), D*~ reconstruction due 
to kinematic fitting and particle identification (5.0%), and Monte Carlo statistics (5.0%). 
Systematic uncertainties which differ for the two modes, and which we quote in parenthesis 
for B^^ D*~ p p 11^ and B^^^ D*~ p n, respectively, are: tracking, 1% per track (6.6%, 
4.6%), proton identification criteria (8%, 4%), and n-body versus two-body for decay kine- 
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FIG. 4. (a) Antiproton reconstruction efficiency. Squares represent Monte Carlo efficiency 
and triangles that of the data, (b) Monte Carlo antineutron momentum spectrum in simulated 
B^ -^ D*~pn decays. 



inatics (5%, 3%). No statistically significant A baryon contribution to the D*^ p p tx^ yield 
was found. We place a systematic uncertainty due to a possible A baryon contribution to 
the i?"^ D*~ p p n^ signal. Also, we assign a 5% uncertainty on the yield of B^^ D*~ 
p n due to the possibility of A baryons with a missing vr distorting the background shape 
in this mode. The systematic uncertainty for antineutron identification is estimated to be 
15%. The quadrature sum of all systematic uncertainties is 15% for B^^ D*^ p p 7r+, and 
19% for 5°^ D*' p n. 

In conclusion, we have made the first observation of decay modes of the B^ of the type 
B -^ DNNX, which may contribute substantially to the total observed B -^ Baryons rate. 
We measure the branching fractions B (5°— > D*~ p p tt"*") = (6.5l};2 ± 1.0) x 10~^, and B{ 



B"^ D* 



p n 



U.dtro ± 2.7) X 10" 
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